Multipartite Entanglement in a One-Dimensional Time Dependent Ising Model. 
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We study multipartite entanglement measures for a one-dimensional Ising chain that is capable 
of showing both integrable and nonintegrable behaviour. This model includes the kicked transverse 
Ising model, which we solve exactly using the Jordan- Wigner transform, as well as nonintegrable 
and mixing regimes. The cluster states arise as a special case and we show that while one measure of 
entanglement is large, another measure can be exponentially small, while symmetrizing these states 
with respect to up and down spins, produces those with large entanglement content uniformly. We 
also calculate exactly some entanglement measures for the nontrivial but integrable case of the kicked 
transverse Ising model. In the nonintegrable case we begin on extensive numerical studies that shows 
that large multipartite entanglement is accompanied by diminishing two-body correlations, and that 
time averaged multipartite entanglement measures can be enhanced in nonintegrable systems. 

PACS numbers; 03.67.Mn,05.45.Mt 



I. INTRODUCTION 



^ . The strictly quantum mechanical property of entanglement has attracted much attention recently, mainly due to 
' its role in quantum protocols such as teleportation, dense coding and other processes that involve transfer of quantum 
^0 information. Entanglement has thus been thought of as a resource for quantum information processing, and perhaps 
Q.^ quantum computing. While there is an understanding of what entanglement is, measures of the same are not so 
obvious, or well established. 

. Entanglement as quantum correlation has also been recently studied with the helpof a slew of well-known models 
' from condensed matter physics, such as the Ising and the Heisenberg models 0, S IS S ■ Mainly, two-body 
, correlations characterized by the concurrence Q have been studied in these systems. Also these were concerned mostly 
^ I with stationary state properties, especially ground states. The entanglement content of a spin-chain, consisting of 
^ I ^ , many spins, could be potentially much more than those that are present in two-body correlations, and nonstationary 
(— I ■ states are of potential interest in small chains, such as those that may be realized in ion trap experiments. The 
' difficulty is in defining proper measures of global entanglement content in such chains. Also it is important to note 
^ . that much of the work has centered around those models that are completely integrable, mostly solvable by the Bethe 
^ '■ Ansatz or by the Jordan- Wigner transform Q . 

^ , Two-body or bipartite entanglement in pure states and its relation to chaos has been investigated more thoroughly 
mainly due to the von Neumann entropy of the reduced density matrices being an unambiguous measure of entan- 
^ . glement. One of the first works to find that chaos leads to larger entanglement production in this case, linked the 
■ classical Lyapunov exponent with the rate of entropy production |lC| . In this case it has been generally found that 
chaos encourages entanglement m 111 mill 111 m, and that complete chaos leads to an universal distribution of 
the eigenvalues of the reduced density matrices giving rise to an universal entanglement that depends only on the 
Hilbert space dimensions |l4l l T5 J . The first study that addressed the role of nonintegrability in many-body entangle- 
ment used the Harper model Il7|. while later works used the quantum baker map |l8| , the Frenkel-Kontorova model 
and disordered spin chains 12011. 
While the relation between entanglement and chaos or nonintegrability is subtle even in bipartite systems, it gets 
even more so in the case of many-body systems. It has been claimed that opposite effects have been observed in 
this case, namely a decrease of entanglement with chaos 0|. However in the case of one-particle states it has been 
observed that the average of all the two-body correlations present in the system does increase with chaos while 
near-neighbour correlations decreases with chaos, where the nearness of the neighbour depends on a kind of quantum 
correlation length Thus it would seem that chaos in these cases can encourage distant entanglement, even of 

a two-body type. However most studies have addressed two-body entanglements, and not global or multipartite 
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entanglements. The exceptions are recent works of Scott and Caves [l8j that make use of a measure due to Meyer 
and Wallach [JJ, called here the Q measure, and indeed show, using the examples of a quantum kicked rotor and the 
quantum baker map, that an increase in chaos entails larger global entanglement. 

Admittedly, global measures of entanglement are only now beginning to be explored and it is likely that the 
various proposed measures quantify different aspects of entanglement in multipartite states, aspects that need further 
elucidation. Wc now briefly recapitulate the definitions of the different entanglement measures used in this paper. 
We emphasize that we will throughout this paper deal exclusively with pure states. 

1. Concurrence: The concurrence in two qubits i and j that are in the joint state pij is given by the following 

procedure 8]: calculate the eigenvalues of the matrix pij plj, where plj = (E) pl^a^ ® , and the complex 
conjugation is done in the standard computational basis. The eigenvalues are positive and when arranged in 
decreasing order if they are {Ai, A2, A3, A4}, the concurrence is C^.j = max (VXi" — ~ \/A3 — ^/Xl^ O). This 
is such that < Cij < 1, with the concurrence vanishing for unentangled states and reaches unity for maximally 
entangled ones. The entanglement of formation of the two qubits is known to be a monotonic function of Cij 
and hence concurrence is itself a good measure of entanglement. In the case of many-qubit pure states, we 
study concurrence between any two qubits by tracing out the others qubits, and studying the resultant density 
matrix. Thus this is a "two-body" correlation. It is known that in typical states multipartite entanglement is 
shared among many qubits rather than in a pairwise manner 18j . We use the following two measures to study 
multipartite entanglement. 

2. Residual tangle and the A^- Tangle: If two qubits i and j are in a pure state the concurrence reduces 

to KV'lo'^ (X) a'^\ip*)\. The state j^/;*) is such that its components in the computational basis are the complex 
conjugates of those of \ip). It was found that the square of the concurrence Tij — Cfj is a more natural measure 
and is now called the tangle We can also define the tangle between one spin (say the fc-th) and the rest of 
the spins if the overall state is pure. This is because in this case, the Schmidt decomposition gives two unique 
eigendirections to the rest of the spins corresponding to those eigenvalues of the reduced density matrix that are 
nonzero. There will be utmost only two such values as the nonzero eigenvalues of the two parts are identical. 
Thus the rest of the qubits can also be effectively thought of as a two-state system. The tangle between spin k 
and the rest, the one-tangle, is t^. (rest) — 4det(pfe), where pk is the reduced density matrix of the fc-th qubit 
(we will also call this simply Tk, not to be confused with the n-tangle introduced below). This was used to define 
a purely three-way entanglement measure in a pure state of three qubits as 



This quantity, called the residual tangle, is independent of the focus qubit, in the above this being the qubit 
numbered 1 ^22,]. and hence stakes its claim as a pure three-way entanglement measure. The construction used 
here to define the residual tangle was generalized to N qubits in Ref. effectively defining a measure of 

multipartite entanglement, the n-tangle as tat — \{iP\ay . This is evidently the tangle for iV = 2, for 

N = 3 this is the residual tangle, while for iV > 3 and odd this vanishes. Thus this measure is used only for 
A'' even, with the exception of A^ = 3. It has been shown to be an entanglement monotone and hence is a 
candidate for measuring multipartite entanglement. It is maximal (unity) for GHZ type states, but can also be 
maximum for states such as the product state of two groups of four spins in the 4-GHZ states. It is of course 
zero for completely unentangled states. It must be noted that the direct generalization of the residual tangle in 
Eq. is conjectured to be positive |23|i and is not the same as the n-tangle for n > 3. 

3. The Meyer and Wallach Q measure: The geometric multipartite entanglement measure Q j21|, has been 
shown to be simply related to one-qubit purities fSJ] , which makes their calculation and interpretation straight- 
forward. It also seems to have the potential for being experimentally measurable. This is defined as: 



From the unit trace of density matrices, it follows immediately that for qubits 1 — Tr(p^) = 2det(pfc). Thus we 
get that 



^1,(23) - n,2 " Ti,3- 



(1) 




(2) 




(3) 



k=l 



This measure is therefore simply the average of the tangle between a given qubit and the rest, averaged over 
this "focus" qubit. The relationship between Q and single spin reduced density matrix purities has lead to a 
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gen eralization of this measure to higher dimensional systems and taking various other bipartite splits of the chain 
|25|. In some ways, for many states, the n-tangle and the Q measure seem to be measuring quite "orthogonal" 
aspects of entanglement, as we see below, eventhough we can and do construct GHZ type of states that maximize 
both these measures. 

It needs to be stressed that it is not true that nonintegrability in general produces more entanglement for arbitrary 
states. There are very simple operators, trivially integrable ones, that can create maximally entangled states out of 
particular unentangled initial states. However, apart from being true only for particular initial states, the entanglement 
will oscillate in time and can be completely destroyed once again. An example is provided by the Hamiltonian 
H = S\® S% of two spin-half particles. With ft = 1 we get, 

exp(-iJS'l ® Slt)\ll) = cos(Jt/4)|ll) - i sin(Jt/4)|00) (4) 

where the states are in the standard 5*^ diagonal basis and |1) is the state with eigenvalue 1/2. The 2-tangle is simple 
to calculate and is 

r=|(Vk>a||^*)P=sin2(Jt/2), (5) 

which follows on substituting the above state. The quantity r is unity for maximally entangled states and is zero for 
unentangled states. Thus after a time t = tt/ J, the spins will be maximally entangled, whereas after twice that time 
they would be totally unentangled once more. Generalizations of such Hamiltonians, and states, to larger spins and 
to larger number of qubits also yield similar results, and is elaborated upon later below. 

In this paper we study issues related to entanglement sharing in spin chains that can range from the integrable to 
the nonintegrable, but which nevertheless involve only nearest neighbour interactions and are translationally invariant. 
This is in contrast to models which have been studied so far, which are essentially single body dynamics, such as the 
Harper or the quantum baker map, that have been mapped onto many qubit systems by means of an isomorphism 
of the Hilbert space. This implies that the interactions need not be nearest neighbour and can in fact involve all to 
all interactions. In particular the model we study is a kicked Ising model of which the kicked transverse Ising model 
is a special case. The kicked transverse Ising model is integrable and we solve it using the Jordan- Wigner transform, 
and thereby study the entanglement generated by this evolution. The zero-field version of this is trivially solvable 
and a class of states that follow in this case have been previously studied as the "cluster states" |2^ . We show that 
while the cluster states have large entanglement as measured by one entanglement measure, Q measure, it has an 
exponentially small (in number of spins) n— tangle. We also show how symmetrizing the states produces those that 
have large entanglement according to both these measures. We also emphasize that the kicking is unlikely to be a 
crucial aspect for the issues discussed here, and on the contrary is more suitable for implementations in say ion-trap 
experiments. 

We numerically study the nonintegrable case and compare it with the integrable one. We find that while nonin- 
tegrability does discourage two-body entanglement, multipartite entanglement is increased on the average. In fact 
this entanglement comes at the cost of two-body correlations. Again, nonintegrability is not required to produce 
maximally entangled states, but produces states that retain large entanglement without disentangling. Thus we study 
the time averaged entanglement measures starting with the vacuum (all spin down) state and study it as a function of 
the strength and tilt of the external field. It is seen that the parameter space corresponding to nonintegrable chains 
is capable of having substantial entanglement. 



II. THE KICKED ISING MODEL 

The model with which we principally study these issues in this paper is a variant |27j | of the transverse Ising model, 
a variant that is at once both dynamically interesting and easier to implement with present day quantum technologies. 
The usual transverse Ising model has been studied in the context of both entanglement and state transport. It is 
an intriguing model that is integrable due to a mapping via the Jordan- Wigner transformation, from interacting 
spins to a collection of noninteracting spinless fermions. The relevance of this model to many physical systems has 
long been appreciated and it is a well studied model, with a quantum phase transition separating ferromagnetic and 
paramagnetic phases at zero temperature as a parameter is varied. The Hamiltonian for L spin 1/2 particles is 

L L 

Hi = Sf^S^^i + S^, (6) 

n— 1 n— 1 

where J is the local exchange coupling strength and B an external transversal field. For J > 25, the system is in a 
ferromagnetic phase with nonzero expectation values of the component of the spin, while for J < 2B the system 
is paramagnetic with vanishing S*^ spin expectation value, the point J = 2B being a quantum critical point. 
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The variant mentioned above involves applying a tilted external field impulsively at regular intervals of time |27l | . 
The operator that evolves states from one application of the field to the next is the quantum map or propagator whose 
spectral properties determine the time evolution. The Hamiltonian is 



oo . L 



H = jY, +B ^ 6ik--)J2 (sin(e)5^ + cos(0)5^) , (7) 



T 

k— — oo 

while the unitary quantum map is (the time T between the kicks sets the time scale and is set to unity): 

U = exp (^-*^E SnS^^+i^ exp (^^^^E + sin(0)5^) j . (8) 

When the field is transverse {9 — tt/2), due to the noncommutativity of the components of the spin operator, the 
above is not equal to exp(— iff/), and gives rise to genuinely different dynamics. However it has been shown that 
this "kicked" transverse Ising model is integrable 27] and there are suggestions to show that it also undergoes a 
quantum phase transition and belongs to the same universality class as the usual transverse Ising model |28| . In this 
integrable model too the key is the Jordan- Wigner transformation, and we solve the problem exactly as opposed to 
the assumption of the thermodynamic limit in Ref . |28j . For < 6 < tt/2 it appears that the model is nonintegrable 
and capable of showing mixing behaviour in the thermodynamic limit [27| . 
Define the following unitary operators: 



UaaiJa) = [] {~^JaS^S^+,) , (9) 

n=l 

L 

UxAB, &) = Y[ exp i-iB (cos(0)5^ + sin(0)5^)) . (10) 

n=l 

Here Lq — L for periodic boundary conditions and is L — 1 for open chains, and 9 is an angle of tilt of the magnetic 
field in the x — z plane. The letter a can be x or z. For the most part we will consider the operator 

U ^U^AJx)U..AB,0). (11) 

This series of unitaries are quantum gates on nearest neighbour pairs of qubits and on individual qubits. Ion-trap 
quantum computing provides one way of implementing the above. The two-qubit operator Uxx maybe implemented as 
phase gates and the single one which involves rotations is implemented via a single Raman pulse. Thus these quantum 
maps maybe experimentally implementable within these architectures in the immediate future. For further details 
and references we refer the reader to 28]. The tilted field changes the character of the dynamics, the Jordan- Wigner 
transformation does not reduce the problem to one of noninteracting fermions and there are features of quantum 
nonintegrability. This has been studied to some extent in the works of Prosen [27| . where he has shown different 
parameter regimes where there is non-ergodic to fully ergodic and mixing dynamics in the thermodynamic limit. This 
model with the tilted field is then one of substantial richness which deserves to be further studied in itself. We will 
use it as a simple and realizable model to study the entanglement issues that were discussed in the introduction. It 
is also worthwhile to mention that time evolution can be done with fast numerical algorithms, with a speed up factor 
of the order of 2^ / L to evolve a state one time step, exactly as the fast Hadamard or the fast Fourier transform. 



III. ENTANGLEMENT IN THE INTEGRABLE CASES 



A. Zero field 



The simplest nontrivial special case of the models in this paper is an extention of what we discussed in the 
introduction to many qubits. Thus we first discuss the set of states: 

\Mt)) ^UMYW"-- (12) 

Here for simplicity we have taken the state with all spins up rather than down. This set of states has been discussed 
earlier and when t/J^ = tt, Stt, 57r, . . . the states are interesting examples of seemingly highly entangled states. 
We say "seemingly" as it is not clear that all measures of multipartite entanglement will be large for these states, for 
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instance we show below that the n-tangle vanishes for these states, when the number of qubits is larger than 3. For 
the case of three spins, such states are local unitarily equivalent to the GHZ state. Expressing the initial state in the 
Sx basis we can easily time evolve and converting back to the standard Sz basis we arrive at: 

I^L(t)) = ^ E expf-^5:(2a,-l)(2a,+i-l)')(g)f^>±i-ll») (13) 

^ a,={Q,l} V fc=l / fe=l \ V2 / 

The states \iPl{t^/ Jx)) are of special interest. For instance for L = 2 we have seen in the introduction that this is 
essentially one of the maximally entangled Bell states. Also up to an overall phase 

\M^/Jx)) = I (|111) - |100) - 1010) - 1001)) , (14) 

which after a local phase change |0) V^1|0), |1) |1), and a tt/4 rotation (Hadamard transform) on each spin 
becomes the GHZ state (|000) + |lll))/\/2. Similarly up to overall phases we get: 

\M^/J.)) = I (10000) - 11111) - llOlO) - 10101)) . (15) 

\M^/Jx)) = J(llllll) - TT (111100) + IIOIOI) + IIOOOO))). (16) 

The operation tt on the states stand for all the five cyclic permutations of this one. For these states we have 
assumed periodic boundary conditions on the spins with Lq = L. We see however that open chains also give rise 
to similarly entangled states. It has been established earlier that these states with L > 3 are not locally convertible 
to generalized GHZ states by means of LOCC. In some sense that has been termed persistence, these states possess 
higher entanglement content than these N-GHZ states or cat states which are jO)®''^ + |1)®^. Persistence is the 
minimum number of local measurements that render the state completely disentangled for all possible outcomes ,26j . 
In terms of a multipartite generalization of the Schmidt numbers, these states seem to again have larger entanglement 
than the GHZ. 

To calculate the Q measure we find the single qubit reduced density matrix pk which is 




Pfc = M 1 (17) 



where the first element is (0|pfc|0) etc., and the angular brackets are expectation values corresponding to the full pure 
state ijj we are interested in. The purity is easily expressed in terms of these expectation values from which we get 
the entanglement measure as 

fe=l fc=l 
For the states under consideration \'ipL{t)) we may explicitly calculate these to get 

(5fe^) = icos^(J.t/2), {S+)^0, (19) 

and hence 

Q{^l) = l-cos\j,t/2). (20) 

Thus this measure of entanglement for this class of states is independent of the length of the chain L, and periodically 
reaches a maximum at t = t^/Jx- Sn/ J^, ■ ■ as indicated earlier, and this maximum is the highest possible. At 
t = 0,2^/ Jx, ■ ■ ■ the state is completely unentangled and therefore in this simple time evolution we have large 
entangling and disentangling oscillations. The periodic boundary condition can be replaced by an open chain, in 
which case the entanglement content as measured by Q is 

Qopen = 1 - COs\jx t/2) - siu^Jx t) , (21) 
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implying again maximal entanglement at times that are odd multiples of ■k/ Jx- Notice that for open chains there is 
marginal dependence of Q on the number of spins, and for L = 2 this simplifies to sin^( J^. t/2), which we have already 
derived as the 2-tangle for this state in the introduction. 

If for these states there is high entanglement content as measured by Q, the two-spin correlations as measured by 
the concurrence is of interest. For L — 2 the (square of the) concurrence coincides with Q, but for higher number 
of spins, we find that while nearest neighbour concurrences persist and oscillate in time, all other concurrences are 
perpetually and strictly zero. Also the times at which the nearest neighbour concurrence vanish are periods when the 
multipartite entanglement content as measure by Q is maximized, indicating that two-body correlations are being 
distributed more globally. To calculate the concurrence between any two spins, at say positions i and j, of the chain, 
we need the two-spin reduced density matrix which is 

{ab\pij\cd) ^ ^ {siS2-..a...b...SL\tp){il;\siS2...c...d...SL), (22) 

Ske{o,i} 

where a and c are fixed states at position i (0 or 1) and similarly b and d are at position j. This matrix can also be 
written in terms of spin-expectation values as 

(00|p.,|00) = ((i - Sf)(i - S^)), (00|py|01) = ((i - 5f )5+), (00|p,,|10) = (5+(i - S])), 
(00|p.,|ll) = {S+S+), (01|p,,|01) = ((i - 5f)(i + 5|)), (01|p,,|10) = (S+Sr), 

(oi|p,,|ii) = {s+i^ + s^)), (io|p,,|io) ^ ((i + sml - s^)), 

(10|P.,|11) = id + S^)S+), = {(l + Sf)(i + S])). (23) 



The rest of the matrix elements follow from Hermiticity of the density matrix. For the class of states given by 
|'0l(^)) can calculate these expectation values in a straightforward manner, exploiting the translational symmetry 
of the states. We get that if j ^ « ± 1 that the density matrix is diagonal, in fact 

Pij ^p^®pJ, j^i±l. (24) 

Here pi and pj are the single spin density matrices as given in Eq. (|17() . Thus there is no concurrence between spins 
that are not nearest neighbours. For the case when j — izLl we get that 

(00|p,,,±i|00) = (01|p,,±i|01) = (10|a,±i|10) = isin2(J,</2), (ll|p,,,±i|ll) = 1 - ^sm^iJxt/2). (25) 
The only nonzero off-diagonal matrix element is 



-I 



(00|a,±i|ll) - — sin(J,t) (26) 

For density matrices such as we have, with all vanishing off-diagonal elements except the corner ones, it is easy to 
find the concurrence in terms of the matrix elements of the density matrix itself. The positive square-roots of the 
eigenvalues of the matrix pi^i±i Pi,i±i, arranged in nonincreasing order are |&| -I- a/ a(l — 3a), a, a, — 16| + y^a(T^-"3a), 
where a = sin^( t/2)/A and b = \ sin( J^; t)/4|. Thus in this case we get, 

a,tti(i) =max(^0, i(|sin(J,t)|-sin'(J,t/2))^ . (27) 

Thus we can explicitly calculate the concurrences at all times, as we see that times at which Q is a maximum namely 
ai t = tt/ Ja;, 27r/ J^, . . ., all the concurrences vanish, including nearest neighbour ones. In fact there is a period of 
time around when Q reaches its maximum that there is no two-body entanglement at all. We get that Ci_i±i = if 
\i&n{Jxt/2)\ > 2 or if t = 2T:k/Jx, k = 0,±1,±2, .... 

In Fig. is shown the entanglement measure Q and the nearest neighbour concurrence as a function of time. This 
figure is independent of the number of spins in the chain, as long as L > 2. The concurrence is dominated by other 
types of entanglement. It has been conjectured that [22| 

n-Y.Ch>^- (28) 
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FIG. 1: The Meyer and Wallach measure of entanglement Q and the nearest neighbour concurrence for the state \^l)4{t)) as 
functions of (scaled) time. Plotted are the numerical (points) and the formula (solid line). Periodic boundary conditions were 
used. 



As we have shown earlier the average of is nothing but the entanglement measure Q, and from translational 
invariance of the states under discussion, this is also equal to any Ti. For this class of states the inequality is easily 
seen to be rigorously true. This difference is interpreted as the generalization of the residual tangle, entanglement 
not present in the form of two-body correlations. In Fig. Q is plotted this residual tangle which is dominated by the 
tangle of individual spins with the others, and not by the concurrence. 

Although both Q and the residual tangle are maximum for states such as / Jx)) the n-tangle measure, as 
previously stated, vanishes for these states. We may calculate explicitly this measure for the states in Ea. l|13|) . and 
we find that 

rjv(V'L) = 1^ sm^{Jx t). (29) 

Thus the n-tangle decreases exponentially with the number of qubits for the cluster state, it seems to be a rare 
entanglement feature, and in particular for the states at t/J^ — tt,2tt,..., such as those in Ea. (|16|l the n-tangle 
vanishes. We note in passing that one class of states for which both the n-tangle and the Q measures are high are 
easily obtained from the states discussed here so far by symmetrizing with respect to the "up" and "down" spins. 
Thus we consider initial states that are N-GHZ states, with the dynamics of nearest neighbour coupling. 

\Mt)) = c/L ^ (|o)^^ + = ^(1 + ®Li<yl\Mt))- (30) 

The last equality follows since the time evolution commutes with the operator ct^ that flips spins in the standard 
basis. For these states (S^) = {S'^) = for all fc, implying that the single spin reduced density matrix is maximally 
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FIG. 2: The residual tangle for the state \ip4{t)) as a function of scaled time. Plotted are the numerical (points) and the formula 
(solid line). Periodic boundary conditions were used. 

mixed, and the measure Q is unity for all time t. The n-tanglc though changes from the maximal value of unity at 
zero time (the N-GHZ states) and oscillates with exact returns to unity at multiples of t:/Jx- 



z^/2sin^/2 



(J,i/2) 



(31) 



Thus for this class of symmetrized cluster states, the n-tangle does not decrease exponentially with the number qubits 
and can have the maximal value at nonzero times. We remind the reader that this measure requires that the number 
of qubits h be even. 
For L = 3 the state 6 is 



1 



2V2 



e*^/2 cos(t/2)) (|000) + |111)) - i sin(T/2) (7r(001 + 110)) 



(32) 



where we have written the scaled time r = 3xt^ that may be simply viewed as a real parameter. For three qubits the 
residual tangle and provides a global entanglement measure [l^. For the state |03(t)), the one and two spin reduced 
density matrices are simply p\ — 1^1^^ where I2 is the two dimensional identity operator and 



/1 + cos2(t) zsin(2r)/2 isin(2T)/2 ~sa?{T) \ 



^12= 4 



sin (r) 



sin2(T) 
sin^ (r) 



V 



-zsin(2T)/2 
~i sin(2r) /2 
1 + cos2(t) / 



(33) 



The other matrix elements of pn follow from Hermiticity of this matrix. Due to translational invariance these are 
the only relevant operators. The spectrum of p\2 is {0,0, 1/2, 1/2}, independent of the parameter r and the spectrum 
of piipii is similarly {0,0,1/4,1/4}. Thus the concurrence vanishes between any two qubits for all values of the 
parameter (time) r. The tangle between one qubit and the other two is T\ = 4det pi = 1, thus the residual tangle is 
Ti — Cf 2 — C'l 3 — 1- Thus we have a continuous one parameter family of states, of which the GHZ state is a special 
case, that have maximal entanglement, as measured by both Q and the residual tangle. Note that in the case of 3 
qubits the residual tangle is also maximized for all time, a feature that generalizes to higher number of qubits, while 
the n-tangle oscillates as indicated above. 
For L = 4 the state is 



\Mr)) 



1 



V2 



cos2(r/2)(|0000) + |1111)) 



2^/2 



sin(T)^(|1100)) 



1 



V2 



sin2(T/2)(|1010) + 10101)). 



(34) 



While (5(04) — 1 for all r, the n-tangle is maximized for r = tt, in which case the state becomes proportional to 
1 1010) -t- |0101), which is local unitarily equivalent to the 4-GHZ state, by say flipping the first and third spins. 
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However for larger number of qubits, the state that maximized the n-tangle is apparently not the N-GHZ state. For 
instance for L — 6, and 8 we get 

106(71-)) = ^ (000000 + 7r(101000+ 100100- 110000) + 1 ^ 0). (35) 
4v 2 

\M^)) = ^ (00000000 + 7r(00010001- 01100110 + 10101010 -00001111 + 01000100) + 1 ^0). (36) 
There are a total of 32 terms in each state and we have temporarily dispensed with the ket notation. 

B. Transverse field 

We now turn on an external field in the transverse direction. This model, the kicked transverse Ising model, has 
been studied recently as noted above and is also an integrable case [23, 13 ' ^'^"^ ^^i*^ Jordan - Wigner transformation 
can be used to diagonalize it. In this case we have 

\4'L{t)) = {UM)UxAB, 7r/2))* l^i(O)) (37) 

where t is an integer time, the number of kicks. We now proceed to diagonalize the operator, indicating the key 
steps. It maybe noted that unlike the treatment in |28j | we do not assume the thermodynamic limit, and in this sense 
the way we solve this problem is also new, though the technique is the same as that for the usual Ising model in a 
transverse field. 

In the kicked transverse Ising spin chain treated here, the Ising interaction is in a;-dircction and the magnetic field 
is switched on at integer times along the z-direction. The first step is to replace the spin variables by Jordan- Wigner 
fermions through a nonlocal transformation 

S+ = exp |] ct c„^ c\ , 5f = Q - i . (38) 

The operators c/ and c\ obey the usual fermion anticommutation rules. The interaction term in Uxx reduces to a 
combination of nearest-neighbour fermion hopping, pair-fermion annihilation and creation terms on a lattice, 

U.. = exp (^-^ {^{c\ - Q)(ct^i + Q+i) - {-ir^cl - cl){A + ci)^ ^ (39) 

where Np = X^i^i is the total number of fermions. The last term is due to the periodic boundary condition. 
The magnetic field term in Ux.z{B, Tr/2) becomes a chemical potential term for the total number of fermions. The 
eigenstates of U will have a definite even or odd fermion number, since Nf commutes with {/, and we can find the 
eigenstates in the two sectors separately. 

Now, the second step is to Fourier transform through, 

exp (i7r/4) , . ^ .,„x 

Cg = "^"^ ' exp {-iql)ci, (40) 

V-t^ 1=1 

where the allowed allowed values for q are (taking L to be even) 

9 = ±|,±|,...,±^^^ iV. even, (41) 

<Z = 0,±^,±^,...,±(^1^,. Np odd. (42) 

The lattice momentum q labels the momentum creation and annihilation operators that also obey the fermion anti- 
commutation rules. The unitary operator U has a direct product structure in terms of these fermion variables: 

U =e-^^l\q>oV^ Np even, (43) 

= e-'^VoV^llg>oV<, Nf odd (44) 
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FIG. 3: The measure Q for the kicked transverse Ising interaction, when the initial state is the vacuum state and L = 10, and 
the parameters are Jx = tt/2, B = tt/3. Shown axe the results of the numerical calculations (points) and using the formula 
(solid line). Periodic boundary conditions are assumed. 



where 



and 



Vg = exp ( -i— 



^-i Y cos(g) (cjc, + cigC-q) + sm{q){cqC-g + dgcl) ^ exp (^-iB{clcg + cl^c-g)) , (45) 



Vo = exp (^-i{B + y)4co^ , V„ = exp (^-i{B - y )4c^^ . 



(46) 



The eigenstates of U are direct products of eigenstates of Vg. The operators Vo and are diagonal in the number 
basis states. For Vg, the four basis states are |0), | ±g) = Cj_q|0), \ —qq>= clgcj|0). The eigenstates of 1^, for g ^ 0, tt 
are given by 



Vg\±q)= e-'(*+^) \±q), Vg\±) = e-'(*+^)e±'^<' |±). 



(47) 



Here the eigenstates |±) are given by |±) = a±{q)\0) + b±{q)\ — qq). Using cos(^g) = cos(B) cos( Jx/2) 
cos(g) sin(B) sin(Jj;/2), we have 



cos(Jx/2)-cos(6'q±B) 
sin q sin B sin( Jx/2) 



(48) 



b fo)=a ( \ ^ s'"(^g) + cos( Ja;/2) sin E - cos g cos B sin( 7^/2) uq\ 

sin(g)sin(J,/2) ' ^ ' 

This then completely solves the kicked transverse Ising model. Let us consider an initial state with m (even) 
fermions |V'(t = 0)) = \li,h---lm) where U denote the sites occupied by fermions (corresponding to Sf, = 1/2 in terms 
of the original spin variables). The off-diagonal matrix element of pi through time evolution with U is 

(5+(i)) ^= (V(t)|e"5:4c„,t|^(i)) =0, (50) 

as the time evolution mixes only states with even number of fermions. The diagonal matrix elements of pi depend on 
{Sf) = {'^{t)\c\ci\'4}{t)) — 1/2. This can be calculated from the time-evolved operator, 

'^S) = cg V' = CgCg- sgn(g) ry, cL„, (51) 
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FIG. 4: (Color Online) The time averaged Q as a function of system parameters for the kicked transverse Ising model. L = 20 
in this case, and the averaging is done over a thousand kicks, by which time the average is stationary. 



where the expansion coefficients are given as 

C, = |a+(g)|^e-"^' + |a_(g)|V*^^ (52) 
77, = a+(g)*6+(g)c-"«'+a_(g)*6_(g)e"^«. (53) 

The diagonal matrix element can be expressed in terms of the Fourier transforms of the above functions, after some 
manipulations, we have 

1 1 " 

(Slit)) = + - ^ |r;,p + ^ \ai - k)\' - HI - k)\'- (54) 

q i—1 

In the above we used two more auxiliary functions defined by 

2 

T]{1) = -^7yqC0s(g/), (55) 

q>0 

^(') = |EC?cos(<zO- (56) 

q>0 

In particular for the initial uncntangled state \iPl{Q)) — |0)®^, as a special case we can calculate {S[{t)) at any site 
using the above. 

(5f (t)) ^ (MowmLm = I E i^'^i' - 1/2. (57) 

Here the q summation extends to both positive and negative allowed values. Hence using translational symmetry the 
entanglement measure Q is given in this case by 

QiMt)) = 4x (1 - x), X = ^ E l^^l' = ^ E l«+(9)a-(9) sin(^, t)f (58) 

As illustrated in the example (Fig. lO) the oscillations of Q are now much more complicated. The advantage of 
having an easily computable formula such as Eq. H58|l is that we can study the entanglement measures as a function 
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of the interaction strength and transverse magnetic fields more comprehensively. In order to do that we time average 
Q over sufficiently long scales and plot this as a function of and B in Fig. This figure shows some interesting 
features, especially the large Q parts which correspond to — tt. Note that both the lines B = and = have 
been discussed previously, the latter case turns off the interaction and produces no entanglement, while the former is 
the zero field case for which the cluster states were realized. 

The case when Jx ~ , B ~ 7r/2 simplifies considerably, as in this case a± {q) = l/\/2 and dq ^ n - q. Thus Q can 
be calculated more explicitly and results in 

QW = {S!HifcV2 fc-O-M,.... (59) 

Thus either the state is maximally entangled by the measure Q or is not at all entangled. As in the zero field case, if 
the initial state is an N-GHZ state itself, according to the Q measure it remains maximally entangled, as in this case 
also {S^) = {S^) = for all times. As in that case the n-tangle measure is now significant, although not maximal in 
general. In the case B = Tr/2, Jx = tt both the Q measure and the n-tangle are unity and represent highly entangled 
states, which appear to be in the nature of cluster states discussed previously for the zero field case. Incidentally, this 
point is also on the critical line Jx = 2B of the (unkicked) transverse Ising model. 

IV. ENTANGLEMENT IN THE NONINTEGRABLE CASE 

We now consider the case when the field is tilted in the x — z plane, that is the unitary operator is a slight 
modification of the transverse Ising case: 

\Mt)) = {Uxx{J.)UxAB,e)f IV-lW). (60) 

The case when 9 is different from both zero and 7r/2, as has been noted earlier, constitutes a nonintegrable model. 
The Jordan- Wigner transformation no longer renders the problem into one of noninteracting fermions. Here we study 
the influence of this the entanglement content of the states |i/'L(t)), again when the initial state is the "vacuum" state. 
Once more the time t takes integer values. Since the Jordan- Wigner transformation does not help, much of the results 
in this section are done purely numerically, with the help of the fast Hadamard transform. 

We start with a given exchange coupling, and strength of the external field, while varying the angle of tilt of this 
field from zero to ninety degrees, both these extremes being integrable. In Figs. (|5I6|I we see the result of this for 
a particular case. We note that the 9 = Q case is integrable and is essentially the zero field case we have discussed 
earlier. In this case the Q measure of entanglement reaches the maximum value of unity and drops back to zero 
periodically. With a non-zero tilt angle we see that the while the maximum drops from unity, the propensity to 
untangle also decreases considerably, thereby providing on the average larger entanglement than for the zero-tilt case. 
Increasing the angle of tilt further decreases the typical value of entanglement produced. The n-tangle measure shows 
more complicated behaviour, with an intermediate angle producing states that have a large n-tangle. 

The increase in the average multipartite entanglement as measured by Q is accompanied by decreasing overall two- 
body correlations, as captured by the pairwise concurrence amongst the qubits. This is illustrated in Fig. Q where 
we show the sum of the two-tangles (or the square of the concurrences), between a given qubit and the rest of them. 
Due to translational symmetry the sum is independent of the focus qubit. This figure shows the rather substantial 
concurrences that are present in the integrable cases (both = and 9 — 7r/2), compared to the nonintegrable 
ones. Thus like the GHZ state that has no two-body correlations, such as the concurrence, these appear to be highly 
entangled states with small concurrences. The entanglement present in the state appears to be predominantly not of 
the two-body type. In fact we noted this previously for the cluster states that when Q was the maximum possible the 
concurrences identically vanished. 

Thus it appears that both the Q and n-tangle measures are sensitive to the nonintegrability of the spin chain and 
from this preliminary data it is plausible that entanglement is enhanced on the average. We have found this to be 
the case for other values of the parameters, not shown here. We can hold the angle fixed and vary the magnitude of 
the external field. Both Figs. (|8I9|I are of this kind. In this case it is seen that small values of the magnetic field are 
enough to prevent the states from completely disentangling. Larger fields also bring down the average along with the 
fluctuations, till for sufficiently large fields the chain seems to reach smoothly an entanglement plateau. The n-tangle 
again shows more complicated behaviour, and can be substantially large in comparison with the integrable cases. 

In order to see the effect of the angle and field strength more comprehensively, we again time average the entan- 
glement measures. This averaging is done over a large number (1000) of kicks such that the average is stationary. 
The results of this are shown in Fig. H1U|) . where it is seen that the Q measure increases sharply with the angle for a 
fixed magnitude B of the field, and then decreases smoothly till the transverse field is reached. The sharp increase 
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FIG. 5: The entanglement measure Q as a function of time, for different tilt angles of the external field. The parameters are 
= 0.1, B = 0.1, L = 10. 
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FIG. 6: The n-tangle as a function of time, the parameters and tilt angles of the external field are the same as in the previous 
figure. 
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FIG. 7: The sum of the two-body tangles as a function of time, for various tilt angles of the external field. The parameters are 
= 0.1, B = 0.1, L = 10. 
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FIG. 8: The entanglement measure Q as a function of time, for various magnitudes of the external field. The parameters are 
=0.1, e = n/4, L = 10. 



is observed in the case B ~ J^, while smoother behavior is seen otherwise. The n-tangle measure shows similar 
characteristic, except that in one case the transverse field case too has a high average entanglement value. 

We next study the time averaged entanglement measures as a function of field strength and tilt, for a fixed exchange 
coupling Jx. The averaging is done over large enough times to ensure stationarity of this quantity, and is shown in 
Fig. (|ll|l . Only six spins are considered here as for each field configuration time evolution is done one thousand times, 
before calculating the average. However the case of larger number of spins is qualitatively similar. The principal 
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FIG. 9: The n-tangle as a function of time, for various magnitudes of the external field. The parameters are same as that of 
the previous figure. 




features seen for the Q measure is that there is enhanced entanglement for both small, nonzero, field strengths and 
tilt angles. The sharp transition at B — Jj. is seen as a fold in the surface plot of this figure. The high entanglement 
spots fall in approximate hyperbolas in the B — 9 space. The time-averaged n-tangle is also shown in Fig. (|ll|l . where 
the hyperbolic region of high entanglement is also visible, but not so close to the small field and tilt angle values as 
for the Q measure. 

From the results presented so far it appears that entanglement can be enhanced in nonintegrable regions of the spin 
chains, but there could be integrable regions such as for zero tilt angle case which could produce large entanglement. 
We have not shown results for the residual tangle in these cases, as this measure is practically identical to Q, this in 
turn implying that the sum of the concurrences is nearly vanishing. In other words two-body correlations as measured 
by the concurrence is a rare commodity in these spins chains. More work needs to be done, especially with different 
initial states, for a better understanding of the implications of nonintegrablity on the entanglement in spin chains. 
The kicked transverse Ising model in a tilted field is a natural example to explore this further. 



V. SUMMARY 



We have studied the Q, the n-tangle, the residual tangle and concurrence measures for a spin chain that is capable 
of showing both integrable and nonintegrable behaviours. The model is the kicked Ising model, kicked with a field 
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FIG. 11: (Color Online) The time averaged Q and n-tangle as a function of external field parameters for the kicked transverse 
Ising model. Jx = 7r/4 and L = 6 in this case. 



that could be transverse or tilted to the exchange coupling direction. The integrable cases correspond to the zero, 
parallel and transverse fields. In the zero or parallel cases the states generated from the vacuum state are essentially 
the "cluster" states, for which we have derived the entanglement measures and shown that while the Q measure is 
large, the n— tangle measure can be exponentially small and the concurrences can vanish. We also point out that 
symmetrization produces highly entangled states that are capable of both large Q and large n-tangles. 

In the case of the transverse field, we solve the time evolution by means of the Jordan- Wigner transformation 
exactly. This enables calculation of many quantities analytically, of which we have displayed the Q measure and 
pointed out the combinations of field strength and exchange couplings that lead to states with large entanglement. 
The Jordan- Wigner transformation docs not help in the case of the tilted field and is an nonintegrable case that has 
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been previously studied from a fidelity point of view. We have studied this case numerically and shown that time 
averaged entanglement can be enhanced in the nonintegrable cases, however it is quite likely that this entanglement 
is not in the form of two-body entanglements. A more detailed study of the nonintegrable case needs to be carried 
out to fully assess the impact of nonitcgrability on multipartite entanglement. 

The entanglement measures Q and the n-tangie have been calculated for random states and it has been shown 
using quantized chaotic maps that these are realized for states evolving under conditions of quantum chaos |lSl l25| | . 
The random state entanglement measure Q for instance is an overestimate for the kicked Ising model even with a 
tilted magnetic field, most likely indicating the effects of translation symmetry, placing strong constraints on the 
"randomness" of these states. Future directions are many, including a more detailed study of states that been shown 
here to have both large Q and n-tangles, especially from an information theoretic viewpoint. Another, is the evaluation 
of the issues studied here with other multipartite entanglement measures, for instance the distance to the nearest 
completely unentangled state p9| . 
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